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Stochastic wake field particle acceleration in Gamma-Ray Bursts 
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ABSTRACT 

Gamma-Ray Burst (GRB) prompt emission can, for specific conditions, be so powerful and 
short-pulsed to strongly influence any surrounding plasma. In this paper, we briefly discuss 
the possibility that a very intense initial burst of radiation produced by GRBs satisfy the 
intensity and temporal conditions to cause stochastic wake-field particle acceleration in a 
surrounding plasma of moderate density. Recent laboratory experiments clearly indicate that 
powerful laser beam pulses of tens of femtosecond duration hitting on target plasmas cause 
efficient particle acceleration and betatron radiation up to tens of Me V . We consider a simple 
but realistic GRB model for which particle wake-field acceleration can first be excited by a 
very strong low-energy precursor, and then be effective in producing the observed prompt X- 
ray and gamma-ray GRB emission. We also briefly discuss some of the consequences of this 
novel GRB emission mechanism. 

Key words: gamma-rays: bursts, plasmas, acceleration of particles, radiation mechanisms: 
non-thermal 
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Figure 1. Oscillation and radiation produced by a relativistic electron oscil- 
lating in an ion channel 



1 INTRODUCTION 

Gamma-Ray Bursts (GRBs) are major explosive phenomena of our 
Universe in need of an explanation. The prompt GRB emission in 
the hypothesis of beamed radiation (with opening angle 9j = 0.1) 
is a very powerful gamma energy radiation with typical emitted en- 
ergy in the range of 10''' erg lasting a few tens of seconds and with a 
peak photon energy Ep of a few hundreds of keV. The compact na- 
ture of the GRB progenitors and the observed photon energy spec- 
tra most likely require both photon propagation transparency and 
relativistic motion of the emitting source. The evolution of a typical 
GRB can be for simplicity divided in two different phases: the first 
phase occurs close to the progenitor and energizes and collimates 
the emitting particles. This first phase may be related to the emis- 



sion of quasi thermal or Comptonized photons (Meszaros & Rees, 
2005). A second step can lead to efficient particle acceleration far 
out from the central source, and eventually to efficient radiation by 
synchrotron emission in internal and/or external shocks. 

In this paper, we investigate a model of GRB production and 
evolution that assumes: (1) the existence of a thermal-like precursor 
of very large intensity following the energy release from relativis- 
tically moving shells or hydrodynamic fronts originating from the 
GRB compact object; (2) the radiative interaction of such a pre- 
cursor (with a given power-spectrum as seen by an observer at rest) 
with a surrounding (initially not-ionized) gas shell at rest of density 
n ~ lO'cm"' at the approximate distance R ~ 10'^ cm; (3) propaga- 
tion of the radiative front within the shell, with possible activation 
of wake-plasma particle acceleration at the outwardly propagating 
radiative shock; (4) wake-field particle acceleration producing a 
generation of non-thermal electrons whose kinetic energy distri- 
bution function extends up to > 100 Me V; (5) betatron radiation 
of these energized electrons producing the observed non-thermal 
prompt GRB high-energy emission in the X-ray/gamma-ray energy 
ranges. The existence of an amount of material corresponding to a 
Thompson optical depth 1 is derived by the observation made by 
BATSE of late bumps in the GRB light curve (Connaughton 2002), 
interpreted as a Compton reprocessed emission from dense circum- 
burst material at /? = 10'^ cm (Barbiellini et al., 2004). The physical 
origin of this material is at present subject of very intense simula- 
tions of the WR wind, where the wind could be compressed by the 
preceding supernova explosion. In this paper we would like to con- 
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sider only the experimental evidence. In this scenario, we assume 
that the energy attributed to the precursor (or to any electromag- 
netic signal entering the plasma) is Eyje where Ey is the energy 
in the prompt emission and e is the conversion efficiency of the 
precursor energy in the prompt emitted radiation. Since the photon 
component of the precursor has quasi-thermal energy, the absorp- 
tion is very large; the highest energy component of the prompt is 
absorbed at a level of 70%. Our GRB emission model does not im- 
ply a relevant role for the intrinsic or MHD-driven comoving mag- 
netic field in the ejecta and/or radiative fronts. From this point of 
view, it provides an emission scenario quite different from those as- 
sumed in current theoretical models of prompt GRB emission such 
as the synchrotron shock model (Tavani 1996; Lloyd & Petrosian 
2002), the jitter radiation model (Medvedev 2000), the diamagnetic 
relativistic pulse accelerator (Liang et al. 2003), and the saturated 
Comptonization model (e.g., Liang et al. 1997). 

The very recent technical progress on high power laser beams 
has provided experimental evidence supporting the concept of 
wake-field particle acceleration in a plasma excited by an elec- 
tromagnetic wave of duration comparable or shorter than the in- 
verse of the plasma frequency (Tajima and Dawson 1979). A laser 
beam with a duration of 30 femtoseconds and a power of 50 Ter- 
aWatts, illuminating a gas jet 1 mm long with a density of 10" 
particles per cm^^ is able to accelerate electrons to an energy of 
few tens of MeV. An X-ray beam of a few keV generated by the 
electron radiation loss in the plasma was observed by Ta Phuoc 
et al. (2005) with a characteristic angular distribution resembling 
that of synchrotron radiation (this motivates the qualification used 
in the literature of "Laser driven synchrotron radiation"). Previous 
connections between the laboratory plasma physics and the astro- 
physical plasmas have been suggested by Chen et al. (2002), which 
proposed a wake field acceleration excited by the Alfven shocks, 
producing Ultra High Energy Cosmic Rays (UHECR). Simulations 
using three dimensional particle-in-cell code, have been recently 
done by Ng & Noble (2005) studying the dynamics of a neutral, 
narrow, electron-positron jet propagating through an unmagnetized 
electron-ion plasma. In this paper, following the general hypothe- 
sis already suggested by Tajima & Dawson (1979), the processes of 
acceleration and radiative losses induced by ultra short laser pulses 
on the gas target, validated by recent experimental data, are sug- 
gested to be the source of the GRB prompt emission. In our case, 
the laser role is played by GRB precursor photons emitted by the y 
and e* shells collimated in an angle Oj and moving with high rela- 
tivistic velocity and the plasma density of the target gaseous nebula 
is assumed to be on the average 10' cm"^^ over the region producing 
the prompt emission at a distance from the central source of about 
10'^ cm. 

In order to activate the coherent wake-field particle accelera- 
tion in the target plasma in the laboratory, where n^ ~ 10" cm"^, 
the threshold power surface density is 

a-w, ~ 10'**Wcm-2 



condition (2). In particular, this is the key point in the analogy 
GRB -Lab plasma acceleration is that the coherent radiation of lab- 
oratory lasers can be translated into non-coherent (stochastic) ra- 
diation from the precursor if a price is payed in terms of power 
density. The critical parameter in the coherent Lab plasma acceler- 
ation is the power surface density threshold. We show below that 
awi scales linearly with particle density. Its value is assumed from 
the experiment: 



o-M/, = 3 X 10'"( 



10", 



— )Wcm- 



(1) 



In the hypothesis that the GRB precursor carries an energy 
fi'p,- = Eyje and assuming ?p, ~ tg^h ~ 30/(1 +z) s the power density 
on the surface R^6- 
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>2x 10'" W cm 
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which exceeds by several order of magnitude (around 10*) the 
power density threshold (scaled at n = 10' cm"'). 

The threshold to initiate stochastically the wakefield mecha- 
nism is equal to the coherent power density multiplied by the factor 
-sjtpr/Tp that is is the probability to produce by chance the coherent 
WFA, independently in time and space. 

In other words, GRB not-coherent prompt emission is much 
more powerful than the coherent lab radiation: we assume that the 
loss of coherence is balanced by the enormous amount of power per 
square centimeter as well as the possibility of stochastic resonance: 
the wake field acceleration threshold is overcome. 

FigureQfrom Ta Phuoc et al. (2005) shows schematically the 
photon emission principle. If the plasma is excited by the precursor 
photons, the excited region and the emission of the radiated photons 
are likely to be in a cone of angular aperture equal to 9j starting at 
a distance of 10'^ cm from the central compact object. 

Stochastic acceleration and deceleration are in principle re- 
lated by plasma effects. However, since they are statistically inde- 
pendent, the two processes can be separated without loss of gener- 
ality. In the following, the acceleration phase is not treated in detail. 
On the basis of the experimental results of Ta Phuoc et al. (2005) 
it is assumed that a power law energy spectrum of electron and 
positron is produced. The accelerated beam enters in the plasma re- 
gion with constant gas density and radiates in the plasma, following 
again the results of Ta Phuoc et al. (2005). 

The value of the plasma density in the GRB environment is 
assumed following Barbiellini et al. (2004). Because of the relation 
among the precursor and the radiating electrons, to reproduce the 
GRB data a stochastic treatment of the deterministic formulae of 
Ta Phuoc et al. (2005) is applied. 

In the laser-driven betatron radiation the emission follows the 
scheme of Ta Phuoc et al. (2005)[fig. 1]. The angular distribution 
of the emitted radiation is dominated by the wiggling angle of the 
electron, 

„ K 



Furthermore, assuming a coherent precursor emission the two 
conditions need to be satisfied: (1) the radiation has to be suddenly 
delivered (in the form of pulses or micro-bursts) to the target gas 
within ultra-short timescales t,, (that for our conditions need to sat- 
isfy the condition r,, < 10"' s in the observer/nebula frame); and 
(2) the radiation intensity delivered within these ultra-short pulses 
has to be sufficiently large to activate microscopic plasma electron- 
ion separation. In the following, we assume that the GRB initial 
fireball producing a very intense radiative precursor is satisfying 
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the oscillating frequency, where aip = yjne^ Keam^) is the plasma 
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fundamental frequency and tq that could be estimated as: 

C _ Pip 

is equivalent to the impact parameter of the electron into the pos- 
itive plasma column (e is the electron's charge). Interpreting Ap = 
clojp as the screen distance of a charge inside the plasma, if the 
charge is moving relativistically Ai, = c/a)/, and r^ represent respec- 
tively the longitudinal and transverse screening distance. 



2 SCALING RELATIONS IN COHERENT PLASMA 
ACCELERATION 

The key parameter in the coherent wakefield acceleration (CWFA) 
is the power surface density ctw needed to produce the complete 
charge separation over a lenght Ap on a timescale r,, = w^'. 

The energy volume density p^ for complete separation could 
be espressed as: 

a^ e'n^Al 
p, = E^eo = — = 

where o" is the charge surface density. The quantity p^ scales lin- 
early with n while the energy surface density cr^ = p^Ap varies as: 

(Tf = oc n ' . 

The key parameter a^ = cr^ojp then scales as the plasma particle 
density n. 

Using the relation among d and the plasma quantities 6 = j- ~ 
1/7, since rg = Ap/ yfy it is possible to derive: 

yrl oc rT^ (3) 

Let us consider an electron radiating with the typical formula: 

dE 



,dy eV I ,2 

— = —m c — = T|a±l 

dt dt 6neQC^ 

where a^ is the transverse electron acceleration 



(4) 



a_L = roOJl = ro 



2y- 



The energy loss by the electron in one oscillating period is: 

dy 4 



dt 



oc y-rlcul 



with A? = (i)j^'. At the equilibrium Ay oc y^^r^n^'^ = const. Recall- 
ing Eq|3| we obtain: 

and, hence, y oc jr^l^ . Furthermore, using again the relation yr^ oc 
rT^ it is also possible to derive that ro scales as «"'''. Another way 
to obtain the previous relation is to consider that the restoring factor 
K does not depend on the density n so that: 

K oc (jjp y''- ro oc n''^ n"''* ro = const 

obtaining that ro oc n"'/^. 

Notice that the maximum electron Lorentz factor y„,ax mea- 
sured by Ta Phuoc et al. (2005) is, roughly y^ax ~ 200. For a den- 
sity of 10' cm"^, the electrons can be accelerated up to: 



rm.x ~ 200 



10^ 
10' 



19 \ 1/3 



In case of low densities, wake field acceleration can be as ef- 
fective as producing high energy cosmic rays and, consequently 
gamma rays. 

Following Ta Phuoc et al. (2005) and using the fundamental 
frequency of radiation 10 f in the case the plasma acts as a wiggler 
we have 

,2 
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We obtain an estimate of the peak energy of the emitted radi- 
ation Epeak: 

Epeak = 'hci)f= -Ticy^ — 

The energy distribution of the accelerated electrons is 

my) = Noy-- 

assuming the power law index again from the experiment (Ta Phuoc 
et al. 2005). 

The peak energy computed as an average on the electron spec- 
trum is then: 

\^peak/ ~ "7 ' TTminymax 

4 Ap 



3 GRB GAMMA EMISSION IN THE STOCHASTIC 
WAKE FIELD ACCELERATION REGIME 

We address here the GRB prompt radiation produced by electrons 
radiating in a plasma excited by a very intense precursor playing 
the role of a "laser". The basic formulae of laser-driven betatron 
radiation are assumed taking into account the stochastic nature of 
the acceleration and the radiation phases. Let us now apply these 
results to the case of a plasma accelerated by the GRB precursor 
photons. We will treat the motion of the relativistic electrons ra- 
diating in the plasma. The fundamental difference with the above 
treatment resides in the stochastic nature of acceleration and radia- 
tion in the case of the GRB emission. 

If the power surface density cty/ at any position of the electrons 
producing the GRB exceeds the value measured by the experiment 
(Ta Phuoc et al. 2005) scaled at the GRB density and multiplied 
by the stochastic factor, the radiated photons of the GRB will be 
emitted from any accelerated electron provided that the stochastic 
angle of the electrons 0eff remain below the jet opening angle 9j. 
The elementary angle of the electron oscillation is: 

e,= '-^ 

A,, 

The stochastic motion of the electrons in the plasma results in an 
effective solid angle spanned by the electrons: 

where we have introduced the adimensional quantity A': 



N ■■ 



R 



/2r^„ 



4.3 X 10" 



where R is the distance travelled by the electron inside the excited 
plasma. This will be used to estimate the stochasticity of the elec- 
tron geometrical path. 

An electron with energy E = ynu? wiggling inside the plasma 
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with density n will contribute to the GRB detectable emission for a 
path Rg if: 



[e,s{Rf,,y,n)f^e] 



(5) 



Putting esplicitly the dependence on the path Rg the previous 
equation could be written as: 

Re^(2y)-'' < 9j 
^P 

This relation expressed in terms of the plasma density and us- 
ing the scaling laws for rg, Ap and y derived previously could be 
written as: 

Ro(y) = 2y/2Ap9]y"^ 

The equation |5| introduces a threshold for the electron emis- 
sion over a distance R. The threshold Lorentz factor 7, is obtained 
requiring that even in the case of a large energy loss, the electron 
trajectory still remains within the energized cone. 

The effective energy loss in the stochastic acceleration regime 
is: 

,0.5 



dy 
dt 



. W7 



2 4"'ft lR\ 



where for simplicity we have introduced the classical electron ra- 
dius r„ 



and where we have multiplied the coherent energy loss by the 
stochastic factor s/R/A. This factor is our best estimation of the 
effect of the turbulence of the circumburst medium. Expressing the 
previous relation in terms of y, and substituting dt = c dR, we ob- 
tain: 



dyy ' = Yir^Ap 



-5/2 ,r,3/2 



dR' 



Jl* 






From this relation we estimate a typical "energy loss" distance R^ 
equivalent to the path where the electron looses half of its energy. 
Integrating the previous equation we get: 

\2/3 



and 



^.~4.1;t„|^| r'/« 



R, - „ . -1/3 l^p 



. ^ 2.8 y 

At \ r. 



The value of 7, is then obtained by requiring 

Reiy,) = ^.(7,) 
obtaining 



7, = 1.2 



2/7 



-6/7 



The opening angle of the jet 6j imposes a threshold on the energy 
of the emitting electron wiggling with an effective angle so that the 
mean value of the peak energy has to be computed over the electron 
spectrum integrated only between 7, and 7max- 



Taking into account the stochastic nature of the phenomenum 
due to the external density turbolence, we can calculate the pre- 
dicted value for <£'peak)eff in the stochastic regime, simply multiply- 
ing the coherent value of (iipcak) by the stochastic factor y/RJA^. 
After an integration between 2y, and 7,, assuming that only these 
electrons contribute to the observed radiation, we obtain: 



1/3 






<''■'■>■• -''O'tM'"'''" 



-8/7 



^peak/eff 

From the previous relations, we have: 

<£'pcak>cff ~ 25 keV — 0. 

Assuming the density n = 10' cm"^^ obtained by the analy- 
sis of the Compton tail phenomenum (Barbiellini et al. 2004), and 
9j=0A (dj x 6°) we obtain in the GRB source frame: 

<£peak>eff « 350 kcV 



4 SPECTRAL - ENERGY CORRELATIONS 

Noticing that the quantity 7, introduces then a dependence of 
(£pcak> on dj, we use the above derived formulae for the peak energy 
of the emission and the effective angle of the radiation to derive an 
estimate for the recently obtained Amati (Amati et al., 2002) and 
Ghirlanda (Nava et al. 2006) relations. 

Considering the electrons with an energy spectrum similar to 
that derived experimentally in the coherence regime, with Lorentz 
factor between 7, and 27,, the mean useful electron energy is: 

^jyy'^dy , 

Ey = Nome— = 2\n2Nomc'y, 

j y-^-dy 

If each electron of energy greater than mc^y, emit ny of mean 
energy (iipeak) then the energy conservations requires: 



mc (7) = ny<£pcak> 



he A 



1/3 



(21n2)mc^7, = n,^l-^ y*" 

The energy conservation requires tiy x y'^ ' oc (£peak) ''"*■ 
The total energy Ey being composed of the emission by the A^o 
electrons results 

NomC'iy) = Ey = Nony(Ep^,,t) = Afo'!y(£peak>''"'(£peak>''"' 

Obtaining: 

Ey=NQX const X (Epe^k)^'* 
If the isotropic energy £1,0 is introduced: 

Inverting the previous equation we obtain iipeak °^ 1^^^ as outlined 
in the so called "Amati" relation. 
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5 CONCLUSION 

We studied the properties of a novel meciianism for prompt GRB 
radiation based on plasma laser- wake acceleration and radiation in- 
duced by an intense precursor irradiating a static gas nebula. The 
stochastic radiative emission of the energized plasma produces a 
photon emission with angular distribution dominated by the elec- 
tron trajectories. Under general conditions applicable to our GRB 
model (precursor, static cloud, laser-wake acceleration conditions) 
we showed that there is a relation between the observed jet prop- 
erty (dj) and an electron energy threshold. In our mechanism the 
energy of the emitting wiggling electrons is the typical energy of 
the laser-wake betatron radiation. Furthermore, we derived a rela- 
tion between source properties such as the coUimation angle and 
the emitted energy. 
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